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ABSTRACT 

The distance scale to cosmic gamma-ray bursts (GRBs) is still uncertain by many orders of magnitude- 
however, one viable scenario places GRBs at cosmological distances, thereby permitting them to be used 

suremenTs of the Hnhh ° g . eX ^ Sloa over a significant range of redshifts z. Also, several recent mea- 
\aJ^i ai?l the Hubble constant H 0 appeanng in the referred literature report values of 70 80 kms -1 

^ ' ( « U 8 h u there 1S s’gmficant debate regarding these measurements, we proceed here under the 

assumpnon that they are evidence of a large value for H 0 . This is done in order to investigate the addi- 

^ ° n C0Sm0 g ) 1C L r d f th3t Ca " be ° btained under this ^POthesis when combined 

W u/ tb ^ age the universe and the brightness distribution of cosmological gamma-ray bursts. 

We show that the range of cosmological models that can be consistent with the GRB brightness dis- 
)" b “ d ° n ’ a Hubble constant of 70-80 km s 1 Mpc ', and a minimum age of the universe of 13 15 Gyr 
onstramed significantly, largely independent of a wide range of assumptions regarding the evolution- 
m ry (h nalare of ^ he burst p ° pulat ' on - Low-density, A > 0 cosmological models with deceleration parameter 
favored 8 1 < 4o < 0 and density parameter a 0 in the range « 0.10 0.25 (Q 0 » 0.2 0.5) are strongly 

Subject headings: cosmology: theory — distance scale — gamma rays: bursts 


1. INTRODUCTION 

Gamma-ray bursts (GRBs) exhibit an angular distribu- 
tion that is consistent with isotropy, yet their brightness 
distribution indicates that the burst sources are either inho- 
mogeneously distributed, embedded in a non-Euclidean 
space, or both (Meegan et al. 1992). Since no known Galac- 
tic objects possess both these properties, and since cosmo- 
logical models explain naturally both the angular isotropy 
and the deviation from Euclidean homogeneity, the hypoth- 
esis that gamma-ray bursts are cosmological (Paczynski 
1991) has gained significant acceptance. 

Several studies have addressed the form of the cosmo- 
logical GRB brightness distribution. Early works (e.g., Mao 
& Paczynski 1992; Fenimore et al. 1993; Wickramasinghe 
et al. 1993) demonstrated that the brightness distribution 
observed with the Burst and Transient Source Experiment 
(BATSE) is consistent with nonevolving, standard-candle 
sources embedded in a spatially flat ( k = 0) universe with 
critical density (D 0 /2 = o 0 = 0-5) and zero cosmological 
constant A, and that these bursts are visible to limiting 
redshifts of order unity. Later analyses (e.g., Emslie & 
Horack 1994) showed that in models in which A = 0, the 
limiting redshift is largely insensitive to the value of q 0 ; 
however larger limiting redshifts are consistent with the 
BATSE data for accelerating ( q 0 < 0) models with nonzero 
A. Models incorporating burst source evolution (e.g., 
Horack, Emslie, & Hartmann 1995; Rutledge, Hui, & 
Lewin 1995; Meszaros & Meszaros 1995) demonstrated 
that higher limiting redshifts are possible for burst popu- 
lations in which the comoving rate density of bursts 
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increases with redshift. Most of the emphasis in these works 
is on “standard” cosmological models in which the cosmo- 
logical constant is assumed to be exactly zero. 

Interest in models with A > 0 has been renewed, 
however, by a recent series of measurements of the Hubble 
constant H 0 in the range 70-80 km s' 1 Mpc -1 , obtained 
through the Hubble Space Telescope ( HST ) and from 
ground-based measurements. For example, Freedman et al. 
(1994), observing the galaxy M100 using HST, report a 
value of H 0 = 80 ± 17 km s -1 Mpc 1 . Pierce et al. (1994) 
obtain a larger value of H 0 = 87 ± 7 km s 1 Mpc -1 
through a study of NGC 4571 utilizing the Canada-France- 
Hawaii Telescope on Mauna Kea. Observations by Tanvir 
et al. (1995) of M96 with HST have resulted in a value of 
69 ± 8 km s 1 Mpc -1 . If correct, when combined with a 
minimum age of the universe of 13-15 billion years, these 
measurements require (e.g., Leonard & Lake 1995) that 
A > 0. 

The value of H 0 in general, and the previous measure- 
ments (as well as others) in particular, is currently the 
subject of debate. It is not certain that H 0 is in the range 
reported by the aforementioned authors, and indeed some 
measurements (e.g., Sandage et al. 1996) yield H 0 values 
significantly less than 70 km s - 1 Mpc - 1 . Our purpose here 
is not to provide a critique or assessment of the validity of 
the many H 0 measurements. 

Instead, we proceed under the assumption that the afore- 
mentioned measurements, as well as others, are indications 
of a large value for H 0 , without significant consideration of 
the arguments for or against either a particular range of 
values for H 0 or a particular published analysis. Under this 
assumption, the combination of a large value for H 0 with 
the independently obtained brightness distribution of 
gamma-ray bursts yields constraints that are unavailable 
from either observation alone. 

We begin with a brief discussion of the mathematics 
necessary to address this problem. This section is followed 
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by analyses of the brightness distribution of the BATSE 3B 
(Meegan et al. 1996) catalog of bursts. Here we obtain x 
probabilities from fits of the observed differential brightness 
distribution N{P) to the distribution predicted by various 
cosmological models. The information obtained from these 
probabilities is then united with the (wholly independent) 
determinations of H 0 and the age of the universe to obtain a 
quantitative estimate of the cosmological model’s ability to 
reliably reproduce data consistent with both sets of obser- 
vations. 


2. ANALYSIS 

2.1. Hubble Constant and Universe Age 
We examine Friedmann dust cosmologies that admit a 
nonzero cosmological constant A. For such models, the age 
of the universe t 0 can be written as 

y(o O' 9o) = T o H 0 

f» (1 +z , )^ 1 dz' 

Jo n / 2 <t 0 ( 1 + z') 3 + (l + Qo — 3 <t 0 K1 +z ') 2 +(ffo-<Zo) 

(1) 

where <r 0 is the density parameter and g 0 is the deceleration 
parameter. For the Einstein de Sitter model (A = 0, o„ = 
q 0 = 0.5), this integral is readily evaluated, yielding y = 3. 

Given a reliable lower limit to r 0 (obtained, for example, 
from detailed stellar evolution calculations), we can use 
equation (1) to obtain the region in the (<r 0 > <Zo) plane that is 
consistent with both the estimated age of the universe t 0 
and current estimates of H 0 . 

To begin the analysis, we must first select an acceptable 
set of H 0 measurements that are commensurate with the 
goal of the investigation, namely, to explore the implica- 
tions of a large Hubble constant in conjunction with the 
brightness distribution of cosmological gamma-ray bursts. 
For the purposes of the analysis, we must assume also that 
these measurements are an indication that H 0 may indeed 
be large and that its value is well estimated by the measure- 
ments in question. To this end, we focus specifically on the 
family of Hubble constant measurements obtained through 
recent HST observations of the Virgo Cluster galaxy M100 
(Freedman et al. 1994; Mould et al. 1995). These are sum- 
marized in Table 1. _j 

The average of the quantities in Table 1 is 78 km s 
Mpc" 1 . We adopt this as our “ best-estimate ” value for H 0 . 
To estimate the uncertainty in this value, we have per- 
formed a computation of the full covariance matrix for 
these seven measurements. This accounts for the common 


TABLE 1 

H 0 Values from M100 Observations 


Value 

(km s' 1 Mpc' *) 


Method 


g0 + 17 Recession velocity 

81 + 11 Recession velocity 

73 + 11 Type II SN distance 

84 + 16 Surface brightness fluctuation 

76 ± 10 Elliptical galaxy velocity dispersion 

82 ± 11 Tully-Fisher relation 

+ 10 Type la SN standard candle 


Reference 

1 

2 

2 

2 

2 

2 

2 


References.— (1) Freedman et al. 1994; (2) Mould et al. 1995. 


sources of error; the distance to M100, found in all 
methods, and the distance to the Virgo Cluster, found in all 
except the Type la SN standard-candle method. 

We obtain a value of H 0 = 78 + 9 km s Mpc . 
Because of the sources of uncertainty common to the 
methods, the combination of a number of different measure- 
ments has only minimal impact on the reduction of the size 
of the overall uncertainty. This value is consistent with 
other recent singular measurements of H 0 , for example, that 
of Whitmore et al. (1995) (78 + 11 km s 1 Mpc" 1 ) from 
observations of M87 in the Coma Cluster, and it is more 
conservative in both magnitude and error estimation tUan 
the Pierce et al. (1994) measurement of 87 ± 7 km s 
Mpc" 1 obtained from observations of NGC 4571 in the 

Virgo Cluster. . , . . , 

A similar set of measurements is required for the age ol 
the universe. We utilize the recent age measurements of 
Bolte & Hogan (1995) (t 0 = 15.8 + 2.1 Gyr), Chaboyer 
(1994) (t 0 = 15.5 ± 4 Gyr), and van den Bergh (1991) (t 0 = 
16.5 + 2 Gyr). These are averaged to determine a best- 
estimate value for t 0 of 15.9 Gyr. The full covariance matrix 
is employed again to determine a measure of the uncer- 
tainty in our best-estimate value of t 0 . However, in addition 
to common uncertainties such as the RR Lyrae distance 
scale, there are others, such as common stellar evolution 
physics incorporated into numerical simulations, that are 
more difficult to quantify. Therefore, to avoid obfuscating 
the main intent of the paper with details of uncertainty 
calculation and estimation, we make the simplifying 
assumption that the off-diagonal terms in the t 0 covariance 
matrix are comparable in size to the typical diagonal terms. 
This assumption, although somewhat arbitrary, is quanti- 
tatively more conservative than the covariance matrix cal- 
culated for the H 0 case, and it yields a result of 15.9 + 2.3 

Gyr. ... 

This estimate of the age of the universe is consistent with 
other recent estimates, for example, the Sommer-Larsen 
(1996) value of 17 + 2 Gyr. Chaboyer (1995) has also esti- 
mated the age of the universe, employing a standard stellar 
evolution model (as well as many modifications) with both 
the Layden et al. (1995) and Walker (1992) RR Lyrae dis- 
tance scales. The Chaboyer (1995) standard values obtained 
with each RR Lyrae distance scale are within 1 o of the 
best-estimate value of t 0 computed above. 

Having obtained best-estimate values for H 0 and r 0 , 
these are then combined to determine a measured value for 
the dimensionless product v’o = t 0 H 0 , resulting in 

y 0 = 1.27 ± 0.23 . (2) 

Figure 1 is a contour plot of y(« t 0 , go) as evaluated from 
equation (1). For reference, spatially flat [(3<r 0 - <Jo — 1) = 
0] models are located along the dashed line labeled k = 0. 
Models with zero cosmological constant are located on the 
dot-dashed line labeled A = 0. The standard Einstein-de 
Sitter model, which possesses both the above properties, is 
located at the intersection of the two lines at which o 0 = 
q 0 = 0.5. The contour for the best-estimate value (2) is also 
explicitly shown in the figure. We observe that the only 
cosmological models consistent with the value Vo ~ 1-27 are 
accelerating (q 0 < 0). This result, which applies for all y > 1, 
is well known. What remains to be determined is whether 
the brightness distribution of gamma-ray bursts observed 
by BATSE is also consistent with any of the cosmological 
models that correspond to the measured value of y 0 . 
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Fro. 1.— Values of y(o 0 , q 0 ) as determined from equation (1). Spatially 
tlat (k = 0) models are located on the dashed line, and A = 0 models are 
located on the dot-dashed line. 


2.2. The GRB Peak Flux Distribution 
We measure the brightness P of a gamma-ray burst in 
terms of its peak flux (photons cm -2 s' 1 ), measured over a 
time interval of 0.256 s in the energy range 50-300 ke V. For 
a given GRB located at redshift z, this brightness can be 
written 


gamma-ray bursts’ luminosity distribution (photons s -1 
keV' 1 ), 


</>(£) = A 0 £ % (5) 

with a > 0, we obtain 


P UI 

4n df 


( 6 ) 


with £ 0 = Iff 4>(E)dE. 

Figure 2 contains both the integral A '(>£) and differen- 
tial N(P) brightness distributions from the third BATSE 
gamma-ray burst catalog (Meegan et al. 1996). The lines 
labeled — 3/2 and —5/2 in the integral and differential plots, 
respectively, indicate the shape of the brightness distribu- 
tion expected from a homogeneous population of bursts in 
Euclidean space. The deviation of the data from this rela- 
tion is apparent, occurring nearly an order of magnitude 
above the instrumental threshold. 

This deviation from Euclidean homogeneity may be due, 
at least in part, to cosmological effects. For a population of 
monoluminous gamma-ray bursts distributed with co- 
moving rate density prescribed by n c (z), the observed rate of 
bursts .4 (>P ) with brightness larger than some value P' 
can be written 


L(zXl + z) 

4 n df 


( 3 ) 


J'(>P') oc 


r 


Mf) 

(1 + z) 


r 2 (z)f(z)dz , 


( 7 ) 


where d, is the luminosity distance, d, = S 0 r(z)(l + z), S 0 is 
the local value of the scale factor, and r(z) is the radial 
coordinate of the burst at redshift z. The functional form of 
r(z) depends on the specific cosmological model in question. 
The quantity L(z ) (photons s' 1 ) is the portion of the bolo- 
metric source luminosity that is accessible to the finite 
energy bandwidth of the detector. It can be written as 

J ^Ezd +2) 

<f>(E')dE' . (4) 

E lU+z) 

Here E 1 and £ 2 are the lower and upper limits to the energy 
visible by the detector (for BATSE 50 and 300 keV, 
respectively). Employing a power-law spectral form for the 


where 


j w ~ j , 

v 2 ff o(l + z) J + (1 + q 0 - 3<t 0 X 1 + z) 2 + (<r 0 - q 0 ) 

( 8 ) 

The factor of (1 + z) in the denominator of equation (7) 
accounts for the fact that one is measuring a rate of burst 
occurrence, and the time between bursts at large redshifts 
will itself be cosmologically dilated by a factor of (1 + z) 
The quantity z(P') is the redshift value that satisfies 
equation 6 for a specified brightness P'. 

The differential brightness distribution N(P ) is defined 
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such that N(P)dP is the rate of observed bursts with bright- 
nesses between P and P + dP m . 


N(P) = 


(9) 


It is straightforward to construct a model differential 
brightness distribution N M (P) from a set of cosmological 
parameters (< 7 0 , q 0 ), spectral index a, burst comoving rate 
density function n c (z), and limiting redshift Zo.s- These hve 
quantities constitute a “ cosmological model.” In the work 
that follows, we shall parameterize the comoving rate 
density as n c (z)oc(l + zf, where 0 is a free parameter A 
value of j8 = 0 corresponds, therefore, to a constant (no 
evolution) comoving rate density. The use of the limiting 
redshift as a parameter in the model assigns a particular 
value of z to a particular peak flux P. This is equivalent to 
attribution of a specific luminosity L (photons s ) to t e 
gamma-ray bursts. In construction of burst models, the 
parameterization of the burst luminosity can be done either 
by assigning a value directly to the bursts (e.g., Fenimore et 
al. 1993), or through the use of z P (e.g., Horack, Mallozzi, & 

Koshut 1996b). . ... », , m . 

The resulting model brightness distribution /V^P) can be 
tested against the observed brightness distribution in 
Figure 2b using the y 1 test. One thereby obtains both a 
value for the reduced yf statistic and a probability P(>X ) 
that a random sample taken from the model would result in 
a distribution of equal or greater disparity than the 
observed data. Low values of P(>* 2 ), therefore, indicate 
that the model is unlikely to be the correct one, as it fails to 
reproduce a model distribution N M (P) similar to the 

observed N(P). , . . 

Figure 3 shows a contour plot of P(>X ) in the ' CT o> 4o) 
plane. In generating this example, we have assumed a con- 
stant comoving rate density of bursts n c (z) and a photon 
spectral index a = 2, and that bursts that are located at a 
redshift of unity produce a peak flux of 0.5 photons cm 
S -1 ( 2() = 1). 

The° compliance of the brightness distribution with all 
decelerating (q 0 > 0) cosmological models noted in the 
introduction is also demonstrated in Figure 3. For this spec- 
tral index, comoving rate density, and limiting redshift, the 
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capability of rejecting models on the basis of P(>x ) 1S 
largely insensitive to the value of both <j 0 and q 0 for deceler- 
ating models. As P(>x 2 ) exceeds 0.5 for all these models, 
the null hypothesis cannot be rejected. 

3. COMBINING THE MEASUREMENTS 
3.1. Probabilities for y = t 0 H 0 
In Figure 1, we presented the computed values ofy(<r 0 , <7o) 
found by evaluating the integral in equation (1) for points in 

the (ff 0 , <?o) plane. We noted also that given the observat ‘ ons 
of Hr, and estimates of the age of the universe, their product 
is measured to be y 0 ± = 1-27 ± 0.23. Assuming a 

(minimal information) Gaussian distribution for the rnea- 
surement distribution variable y', it is straightforward to 
compute for each point in the (<j 0 , <Jo) pl ane a probability 
value 


^(|y' -y 0 l > ly-kol) 


2 P® 

J|y- 


-y'2/2 


sf^Tl J|y-yol/Oy 

I y - ko \ \ 


dy' 


= erfc 


(10) 


Fig. 3.— Contours of P(>x 2 ) ' n the (»o> to) plane for a cosmological 
model with constant comoving rate density n c (z), power-law burst spectral 
index a = 2, and limiting redshift of z 0 . 5 = 1. These probabilities are 
obtained by application of the r test between the differential brightness 
distribution predicted by the model and the distribution obtained by 
BATSE (Fig. 2f>). 


This value represents the probability that a measurement of 
the value of y' will be realized to be equally or more discrep- 
ant from the value y 0 than the calculated value y(a 0 , q 0 ) that 
the model predicts. Points in the (<r 0 , q 0 ) plane that yield 
very low probabilities 0* are unlikely to represent reality, as 
they predict a value y that is inconsistent with the measure- 
ment y 0 at a high level. 

For example, we have seen that the integral in equation 
(1), when evaluated with a 0 = Qo = °- 5 > is e 9 ual j 0 . 2 ^ 3 ’ 
which is approximately 2.6 standard errors removed from 
the measured value of 1.27. Thus, the value of 9 assigned to 
this point in the plane is erfc (2.6/2 ,/2 ) « 9.32 x 10 That 
the inconsistency between the measured value y 0 and the 
value predicted by the Einstein de Sitter model is on the 
order of 3 a has been noted in previous works (e.g., van den 
Bergh 1992), and it is one measure of the relative merit of 
both the value y 0 used here and its computed uncertainty. 

3.2. Combining with the GRB yf Probabilities 
To this point, we have determined the goodness of fit 
between the BATSE gamma-ray burst brightness distribu- 
tion and the distribution predicted by the cosmological 
model in question, as well as the goodness of fit between the 
measured value y 0 and the value y(p 0 , q 0 ) predicted by the 
model. We would like to combine these two pieces of infor- 
mation in a way that treats both of them with equal impor- 
tance and allows for the rejection of a particular 
cosmological model on the basis of its level of disagreement 
with either the brightness distribution or the measured 
value of y. The situation is analogous to the consideration 
of a Galactic geometric burst model that predicts both a 
brightness distribution (e.g., Hakkila et al. 1995) and values 
for the two angular distribution moments (e.g., Briggs et a . 
1996) <cos 0> and <sin 2 b>. 

By adopting the reasonable assumption that the mea- 
surements of r 0 H 0 obtained in part from observations 
carried out with the Hubble Space Telescope are indepen- 
dent of the measured brightness distribution of gamma-ray 
bursts obtained with BATSE on the Compton Gamma Ray 
Observatory, a joint probability Pj can be obtained for each 
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Density r: c r a meter ct 0 

Fig. 4.— Joint probability contours P } for the cosmological model used 
in Fig. 3. These data express the probability that a random sampling from 
the cosmological model would return both a x 2 value more discrepant than 
the one obtained by comparing the model N M (P ) to the BATSE distribu- 
tion, and that a measurement of y would be obtained farther from y 0 = 
1.27 ± 0.19 than the value predicted by the model through equation (1). 
For this set of parameters, many high-density and most decelerating 
models can be ruled out with a high degree of confidence, as they predict 
values of y and/or a brightness distribution that are significantly in dis- 
agreement with observational data. 


point in the (<x 0 , q 0 ) plane by simply multiplying P{>x 2 ) and 
i.e., 


Pj = P(>X 2 ) erfc 


It -Tol \ 

a yo \/2 / 


(ID 


The quantity Pj is the probability that a random sam- 
pling from the cosmological model in question would 
return both a x 2 value more discrepant than the one 
obtained comparing the model to the BATSE brightness 
distribution and a measurement value y' that is farther from 
y 0 = 1.27 + 0.23 than they calculated from equation (1). 

We employ Pj as our statistic to assess a model’s capabil- 
ity to reproduce observations of both the brightness dis- 
tribution and value of y. As such, it is more demanding than 
either of the individual estimators [P( > y 2 ) or alone. It is 
also a significantly stronger test for rejection than a global 
X 2 , computed as the sum of the yf from the brightness dis- 
tribution and the value (y - y 0 ) 2 /a 2 , and it treats the mea- 
surement of y on equal footing with only one bin of the 
brightness distribution rather than the overall goodness of 
fit. 

We stress that Pj is used here only to reject models. 
Models which yield very low values of Pj are unlikely to be 
the actual cosmological configuration found in nature, as 
they predict a brightness distribution and/or value of y that 
is significantly different from those observed. Conversely, 
models yielding large values of Pj cannot be rejected, as 
they reproduce both a brightness distribution and value of y 
consistent with observation. The value of Pj that one 
chooses as the threshold for rejection of a particular model 
is, of course, ultimately a scientific judgment. Figure 4 
shows one example of P/<r 0 , q 0 ) for the same cosmological 
model employed in Figure 3. 


4. RESULTS OF COMBINING MEASUREMENTS 

Figures 5a-5d show the results of combining these two 
probabilities for a variety of assumed cosmological models. 
In Figure 5 a we show the contour levels in the (a 0 , q 0 ) plane 
for which Pj = 0.32, employing models with a nonevolving 
burst population, photon spectral index a = 1.0, and 


various limiting redshifts z 0 ? . Each assumed limiting red- 
shift produces a different region in the plane in which P = 
0.32. The contour labels (1 through 4) denote the values of 
these limiting redshifts. 

We have also indicated with a dashed line the positions of 
spatially flat (k = 0) cosmological models, which have 
(3<x o — do ~ 1) = 0. Models with zero cosmological con- 
stant A are shown with the dot-dashed line. To illustrate 
better the regions of constant P p we have also truncated the 
x-axis at a value of <r 0 = f2 0 /2 = 0.5, the critical density for 
A = 0. Figure 5b presents the same information for cosmo- 
logical models in which a = 2. 

In Figures 5c and 5 d, the value of Pj is lowered to 0.01 for 
a = 1 and a = 2, respectively. As expected, the contours of 
Pj = 0.01 are significantly larger than those for Pj = 0.32. 

For a given point in the ( a 0 , q 0 ) plane, we observe that 
higher limiting redshifts are required to obtain a particular 
value Pj for the a = 1 spectral form than are required for 
a = 2. As an example, the curve Pj = 0.32 found for a = 1 
and a limiting redshift of 2 (Fig. 5a) is similar to that found 
with the a = 2 for a limiting redshift of 1 (Fig. 5b). Physi- 
cally, because the a = 2 spectral form is steeper, its intensity 
in a given energy window is more significantly affected by a 
given amount of cosmological redshifting. Therefore, the 
cosmological bending observed in the integral BATSE 
brightness distribution .4 '( > P) is achieved at lower red- 
shifts, and consequently the limiting redshift for bursts with 
a = 2 is less than that required for a = 1 . 

Figures 6a 6d present similar information to Figure 5, 
but for an evolutionary comoving rate density function 
n A z ) ~ (1 + z). The general forms of the figures are similar; 
however, there is one significant difference, namely, that a 
given value of Pj is found in the nonevolving case with 
lower redshifts than in the evolving case. For example, in 
Figure 6a, no limiting redshift below z 0 5 = 4 results in a 
value of Pj > 0.32, whereas for Figure 5a, values of Pj larger 
than 0.32 can be easily found for these redshifts. Physically, 
this behavior is readily understood. As stronger (larger /?) 
evolutionary forms are utilized, the brightness distribution 
deviates less rapidly from a Euclidean slope with increasing 
redshift, so that higher limiting redshifts are necessary to 
create a distribution consistent with the concave downward 
[in the log . ( (> P) log P plane) curve observed with 
BATSE. 

Our most important result transcends all these dis- 
cussions. By inspection of Figures 5 and 6, we observe the 
Pj contours to be strongly peaked in the lower left region of 
the (<r 0 , q 0 ) plane. In addition, almost any cosmological 
model with a density in excess of the critical density (<r 0 = 
fi 0 /2 > 0.5) will result in a low joint probability P } ( ~0.01 
or less), regardless of the assumption about the burst popu- 
lation (evolution, spectral index, limiting redshift, etc.). 
Exceptions are those rapidly accelerating ( q 0 <$ 0) models 
with steep burst spectra (a ~ 2), and limiting redshifts in 
excess of z 0 5 > 2. The analysis is not as constraining on q 0 
because for a given value of <x 0 , a particular value of Pj often 
occurs for both positive and negative (accelerating and 
decelerating) values of q 0 . Thus, q 0 can be said to be 
“ restricted ” to a given range, but the dynamics of models 
belonging to this range can be quite varied. 

Therefore, we conclude that the combined measurements 
of H 0 , t 0 , and the BATSE brightness distribution make 
cosmological models with the density parameter o 0 near 0.5 
(^o ~ I) or larger quantifiably less likely than low-density 
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Fig 5 -Contours of constant joint probability for nonevolving burst populations. The x-axis has been truncated at the critical density <r 0 -JV2- 0.5 
for ckritv^nthe contoursTu) Curies for P, = 0.32 with a spectral index of a = 1. The limiting redshifts used to obtain a given c^toy a e tebekd on the 
curves, (ft) Similar to (a), but the spectral index used is « = 2. (c) The value of P, is lowered to 0.01 for the same sets of parameters used in (a), id) P, - 

« = 2 . 


models based on their low Pj values, regardless of the value 
of q 0 . The observed value of y 0 and/or the BATSE bright- 
ness distribution N(P) are considerably more disparate 
from those that might be expected from a random sampling 
of these theoretical models. 

5. DISCUSSION OF UNCERTAINTIES AND SYSTEMATICS 

It is unlikely that the peak photon luminosities of bursts 
are the same for each event everywhere in the universe; 
however, several analyses (e.g., Emslie & Horack 1994, 
Hakkila et al. 1996) have shown that the range of lumi- 
nosity for a significant fraction of detected cosmological 
bursts is indeed quite narrow (~ 10), validating the use of 
standard candles as a first-order assumption. 

The use of a power-law spectral form for bursts is also an 
approximation, as burst spectra are known to exhibit evolu- 
tion within a given burst and display curvature (Band et al. 
1993). However, Mallozzi et al. (1996) and Horack et al. 
(1996a) have shown that despite their simplicity, power-law 
spectral forms employed in analyses of the cosmological 
burst brightness distribution yield quantitatively consistent 
answers to those obtained with more complex curved spec- 
tral forms. Further, the values of a = 1 and a = 2 are 
approximate lower and upper bounds to the spectra 
observed for bursts in the range 50—300 keV, and their use 
here is intended to bracket approximately the behavior over 
the range of consistent spectral models. 

One may also use any number of measurements for the 
age of the universe or the Hubble constant. As this is an 


investigation into the implications of a high value for H 0 , 
we must necessarily use measurements in the range of 70 80 
km s" 1 Mpc ‘, as they are commensurate with the goals of 
the paper. However, even within this range there are a wide 
variety of results that could possibly be utilized, each with 
their own advantages and drawbacks. We have chosen to 
focus on recent estimates of the age of the universe in con- 
junction with recent measurements of H 0 obtained from 
HST examination of the Virgo Cluster galaxy M100. We 
have explored the effects of utilizing different H 0 measure- 
ments than those cited here. The use of any other particular 
compatible H 0 measurement in place of those employed 
here will obviously change the computed value of y in 
equation (2) slightly; however, this does not alter substan- 
tially the general conclusions of the paper. 

If we employ a value for H 0 that is at the low end of our 
set of values, for example, the H 0 = 69 km s 1 Mpc 1 mea- 
surement ofTanvir etal. (1995) and a value t 0 = 13.5 Gyr(a 
lower bound to the measurement of Bolte & Hogan 1995), 
one still obtains a value of y = Tq H o ~ !■ This forces strong 
consideration of A > 0 models and virtually eliminates flat 
models with a high density parameter, including the 
Einstein de Sitter model. 

Depending on the assumptions one makes, the shapes of 
the various contours of confidence presented in Figures 4, 5, 
and 6 will also change. Their shapes and level values depend 
strongly on the magnitude of the uncertainty for y 0 . Indeed, 
this was one impetus for utilization of the complete covari- 
ance matrix computation earlier. Future efforts to reduce 
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Fig. 6. — Similar to Fig. 5, but these data are produced using an evolutionary form n f (z) ~ (1 + z) for the burst comoving rate density, (a) No models (<7 0 , 
q 0 ) employing a = 1 result in a P ; greater than or equal to 0.32, regardless of the limiting redshifts investigated. ( b ) Contours of Pj = 0.32, for a spectral index 
of a = 2. (c, d) The value of P ; is lowered to 0.01, and the spectral indices used in these evolving models are as shown. 


the error on the measured value of H 0 to ~ 10% via the 
HST key project (see Mould et al. 1995) will further restrict 
nonrejectable models in the (<r 0 , q 0 ) plane. 

We have noted that the brightness distribution of bursts 
is not sufficiently sensitive to differentiate between A — 0 
models, or models in which q 0 is positive (decelerating). An 
example of this feature was noted in the discussion of Figure 
3, where P(>x 2 ) is in excess of 0.5 for each decelerating 
model in question. Since for many “standard” cosmologies 
the problem is overdetermined, very large values of P(>x 2 ) 
are not uncommon. 

This is one reason for the use of the joint probability Pj 
rather than a different test such as the global x 2 - Indeed, 
when the reduced x 2 value for the model N(P) is signifi- 
cantly less than unity, as is the case for many of the afore- 
mentioned standard cosmologies, the further incorporation 
of a significant deviation in y, treated on equal footing with 
each individual bin of the N(P) distribution, is not enough 
to raise the global x 2 value significantly to a point at which 
the model would be rejected on the basis of this statistical 
test. Therefore, we have chosen the more powerful test Pj 
for rejection of models in the test of the null hypothesis. 

That the BATSE N(P) distribution and values of y can be 
combined in this manner to significantly restrict allowable 
values of (<r 0 , q 0 ) is somewhat fortuitous. If y were measured 
to be significantly less than unity, the brightness distribu- 
tion could not be used to further distinguish among possible 
cosmological models. However, with y in excess of unity, we 
are in a region in which small changes in a 0 or q 0 require 
changes in other parameters to obtain model consistency 


(e.g., Emslie & Horack 1994). Thus, there is significant dis- 
criminatory power in this region, as the low-density, mildly 
accelerating models that cannot be distinguished from the 
high-density, rapidly accelerating models on the basis of y 
produce markedly different N(P) distributions for a given a, 
z 0 5 , and comoving rate density n c (z). 

6. FURTHER ANALYSIS AND DISCUSSION 

Because of the level of uncertainty in the data input to 
this analysis, it would likely not be productive to carry this 
analysis from one of hypothesis testing into one of param- 
eter estimation in any great detail. However, examination of 
those cosmological models which reproduce best both the 
observed value of y 0 and the N(P) distribution does offer an 
interesting revelation. 

For a given rate density function n c {z), burst spectral 
index a, and limiting redshift z 0 5 , the point in the (<r 0 , q 0 ) 
plane that yields the largest joint probability value Pj can 
be considered the “best-fit” cosmological model for these 
input parameters. In Tables 2 and 3, we present the values 
of the maximum probability Fornax) for various assumed 
spectral indices and limiting redshifts, as well as the values 
of the density parameter <r 0 (or Q 0 /2) and deceleration 
parameter q 0 at which these maximum probabilities are 
found. Table 2 displays data assuming a nonevolving 
cosmological rate density function for gamma-ray bursts 
[i.e., n c (z) = constant], and Table 3 employs an evolution- 
ary rate density of the form n c (z) oc ( 1 + z). 

We observe that for <r 0 ^ 0 (i.e., a model universe that is 
not totally empty), the maximum values of Pj are generally 
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TABLE 2 


Values for n f (z) = Constant 


a 

Z 0.5 

Pj ( max ) 

^0 

% 

1.0 

1.0 

0.0073 

0.00 

0.08 


1.5 

0.3096 

0.00 

- 0.52 


2.0 

0.6968 

0.00 

- 0.52 


2.5 

0.7917 

0.04 

- 0.80 


3.0 

0.7906 

0.07 

- 0.92 


3.5 

0.7894 

0.09 

-LOO 


4.0 

0.7411 

0.11 

- 1.04 

1.5 

1.0 

0.1900 

0.00 

- 0.52 


1.5 

0.6918 

0.00 

- 0.52 


2.0 

0.7758 

0.08 

- 0.96 


2.5 

0.7880 

0.15 

- 1.16 


3.0 

0.7491 

0.16 

- 1.20 


3.5 

0.7721 

0.18 

- 1.24 


4.0 

0.6656 

0.18 

- 1.24 

2.0 

1.0 

0.5944 

0.00 

- 0.52 


1.5 

0.7500 

0.11 

- 1.04 


2.0 

0.8326 

0.18 

- 1.24 


2.5 

0.8058 

0.21 

- 1.32 


3.0 

0.8341 

0.24 

- 1.36 


3.5 

0.7846 

0.24 

- 1.36 


4.0 

0.7828 

0.24 

- 1.36 


found for o 0 between 0.10 and 0.25, regardless of the 
assumed spectral index, evolutionary mode, or limiting red- 
shift. These correspond to values of Q 0 between 0.2 and 0.5, 
significantly less than the critical density Q 0 = 1. Each 
model also possesses a positive cosmological constant A 
that can be calculated as (e.g., McVittie 1965) 

A = ( ff o - </o) • (* 2 ) 

Furthermore, each of these best-fit models lies in a region of 
parameter space in which the quantity 3cr 0 — q 0 — 1 (the 
sign of which determines the spatial curvature k of the 
model) is small, in many cases very near zero. 

We encourage readers not to overinterpret these data in 
light of the nature of the assumptions outside the present 
analysis that have to be made. For this very reason, we have 


TABLE 3 


Values for n c { z ) oc (1 + z) 


a 

2 0.5 

Pj ( max ) 

< 7o 

4o 

1.0 

1.0 

1.33 x 10~ 6 

0.00 

1.48 


1.5 

0.0005 

0.00 

0.56 


2.0 

0.0117 

0.00 

- 0.24 


2.5 

0.0668 

0.00 

- 0.36 


3.0 

0.1774 

0.00 

- 0.52 


3.5 

0.3148 

0.00 

- 0.48 


4.0 

0.4342 

0.00 

- 0.52 

1.5 

1.0 

0.0010 

0.00 

0.08 


1.5 

0.0694 

0.00 

- 0.28 


2.0 

0.3302 

0.00 

- 0.48 


2.5 

0.5569 

0.00 

- 0.52 


3.0 

0.6779 

0.04 

- 0.80 


3.5 

0.7869 

0.06 

- 0.88 


4.0 

0.7862 

0.07 

- 0.92 

2.0 

1.0 

0.0488 

0.00 

- 0.36 


1.5 

0.4594 

0.00 

- 0.52 


2.0 

0.6585 

0.05 

- 0.84 


2.5 

0.7894 

0.11 

- 1.04 


3.0 

0.8035 

0.17 

- 1.20 


3.5 

0.7821 

0.18 

- 1.24 


4.0 

0.7094 

0.20 

- 1.28 


not taken the analysis here an additional step further to 
perform detailed calculations of the best-fit values and 
uncertainties for each of the parameters a, z 0 5 , evolutionary 
power-law index, cr 0 , and q 0 . 

However, regardless of the detail of the assumptions, one 
result is clear. If we live in a universe that is older than 
~ 13.5 Gyr and possesses a Hubble constant in excess of 
~70 km s 1 Mpc~\ the cosmological models that best 
reproduce both the brightness distribution of gamma-ray 
bursts and are most compliant with these aforementioned 
numbers possess low-mass density, with £2 0 in the range 
0.2-0.5, and a positive cosmological constant. 

Furthermore, it is interesting to note that for most of 
these models the quantity | 3 <j 0 — q 0 — 1 1 is small, nearly 
zero in some cases. The spatial curvature of a cosmological 
model (k = — 1, 0, -hi) is determined by the sign of 
(3cr 0 — q 0 — 1). Specifically, one has the relationship (e.g., 
McVittie 1965) 



For a fixed H 0> we observe that as the absolute value of 
3a 0 — q 0 — 1 becomes very small, the value of the scale 
factor S 0 becomes very large. Hence, the radial coordinates 
r we can observe are always 1 in these models, even for 
objects that might be seen at exceptionally large redshifts. 
In these models, then, objects cannot be observed at dis- 
tances that are large enough r 1 such that the extent of 
the spatial curvature is appreciable. We refer to these as 
models with “small curvature/’ inasmuch as the “radius” 
of possible curvature is much greater than the extent of the 
observable universe. 

On the basis of other gamma-ray burst evidence not 
incorporated quantitatively into Tables 2 and 3, there may 
be reason to consider further those models with large limit- 
ing redshifts as less likely than those with moderate values 
of z 0 5 . These larger limiting redshifts generally result in the 
“best-fit” models occurring for larger values of <x 0 (and 
therefore n o ). The existence of a large limiting redshift 
requires (e.g., Horack et al. 1995) the — 3/2 power law 
observed at the bright end of the log A r (>P) versus log P 
curve to be a consequence of evolution and/or cosmological 
effects conspiring to produce a brightness distribution that 
looks at the bright end like a homogeneous population of 
sources distributed throughout Euclidean space, although 
in these instances those bursts would be neither homoge- 
neous nor Euclidean. The existence of the —3/2 slope in the 
data is strengthened by the combination of BATSE data 
with observations of Pioneer Venus Orbiter (PVO) 
(Fenimore et al. 1993). In addition, searches for cosmo- 
logically induced time dilation (e.g., Norris et al. 1995) and 
spectral shifting (Mallozzi et al. 1996) have revealed no evi- 
dence of a signal strong enough to be attributed to bursts at 
redshifts of z = 3 or 4, finding instead that the data are 
consistent with limiting redshifts of the order z ~ 2 or less. 
Thus, on the basis of data not incorporated here, there may 
be further evidence to lessen the likelihood of the larger Q 0 
models shown in Tables 2 and 3. 

A significant body of other observations (see, for example, 
Coles & Ellis 1994) indicates, that Q 0 is likely to be in a 
range consistent with that found here. Among them, esti- 
mated abundances of 7 Li, 4 He, 3 He, and deuterium in the 
primordial universe agree with nucleosynthetic predictions 
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(Walker et al. 1991) only if the baryonic component of the 
mass density is in the range 

100 < Q b #o < 150 . (14) 

More recent estimates utilizing big bang nucleosynthesis 
and observed elemental abundances have constrained Q b to 
be less than 0.1 (Copi, Schramm, & Turner 1995). Only the 
postulation of a significant amount of nonbaryonic “ dark 
matter” can result in a large value of Q 0 given these small 
values ofO b . 

Ostriker & Steinhardt (1995) also derive a range of Q 0 
consistent with that presented here through a combination 
of several independent observations of H 0 and the anisot- 
ropy of the cosmic microwave background. 

Inflationary models of the universe (e.g., Guth 1981), 
popular on theoretical grounds, strongly favor zero spatial 
curvature (3 a 0 — q 0 — 1 = 0). The inflationary expansion 
produces enough energy to drive the universe toward a 
critical total energy density 2 <r 0 + (<7 0 — q 0 ) ~ 1, which 
includes contributions from mass density 2<t 0 ( = D 0 ) and a 
cosmological constant A (~<7 0 — q 0 )• Under the assump- 
tion of an appreciable value of A obtained from the mea- 
surements of H 0 , the further strong preference for a 
spatially flat universe mandates a low value for <r 0 = Q 0 /2. 

Given that it is by no means a certainty that gamma-ray 
bursts are cosmological, the concurrence of any cosmo- 
logical model with both the brightness distribution and 
values of y 0 may be treated with some degree of curiosity. 
On the basis of gamma-ray burst data alone, there is no 
reason to favor low-density, A # 0 models over high- 
density ones. Yet it is the N(P) distribution of burst bright- 
nesses which prevents the inclusion of rapidly accelerating, 
high-density models that could otherwise be invoked to 
explain the large value of y 0 . Thus, the combination of the 
BATSE gamma-ray burst N(P) distribution with the mea- 
sured value of y has led to a restriction of the (<r 0 , q 0 ) 
parameter space that is unobtainable from the measure- 
ments of H 0 and t 0 alone. 

7. CONCLUSIONS 

We have presented the first meaningful constraints on the 
cosmological parameters cr 0 and q 0 to be derived, in part, 
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from analyses of the BATSE gamma-ray burst brightness 
distribution. We find that in the context of a large Hubble 
constant (70-80 km s" 1 Mpc - : ) and an age of the universe 
near 15 Gyr, cosmological models that possess a density 
near <r 0 = Q 0 /2 = 0.5 or larger are quantifiably less likely 
representations of reality, as they fail to produce both 
values of y 0 = t 0 H 0 and a burst brightness distribution 
N(P) that agree with the observational data. This conclu- 
sion is independent of a wide range of assumptions regard- 
ing the burst evolutionary rate density, spectral index, or 
limiting redshift. 

Furthermore, the peaked portions of the Pj contours, i.e., 
those models that reproduce most accurately both the 
observed t 0 H 0 value and BATSE brightness distribution, 
lie in a region of parameter space in which <x 0 = Q 0 /2 « 
0.10-0.25, significantly less than the A = 0 critical density. 
Members of this family of cosmological models are all low- 
density, accelerating models with a positive cosmological 
constant A, and minimal spatial curvature ( 1 3cr 0 — q 0 — 1 1 
near zero). 

Coles & Ellis (1994) summarize the observational evi- 
dence favoring a low-density universe, and Ostriker & 
Steinhardt (1995) have argued further on the basis of other 
data, that not only is Q 0 within the range found here, but 
that the universe also is likely to have minimal spatial cur- 
vature and a positive cosmological constant. If recent mea- 
surements of the Hubble constant are substantially correct 
and our estimates of the age of the universe are fairly accu- 
rate, we find that these cosmological models also best 
explain jointly the observed value of y 0 and the BATSE 
brightness distribution. 
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